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Abstract 

Several general constraints are suggested to analyze the pionic valence-state 
wave function. It is found that the present model wave functions used in light- 
cone formalism of perturbative quantum chromo dynamics have failed these re- 
quirements to fit the pionic formfactor data and the reasonable valence-state 
structure function which does not exceed the pionic structure function data for 
x — > 1 simultaneously. Furthermore, it is pointed out that there is a possibility 
to find model wave functions which can satisfy all of general constraints. Also 
we show that there are two higher helicity (Ai + A2 = ±1) components in the 
light-cone wave function for the pion as a natural consequence from the Melosh 
rotation and it is speculated that these components should be incorporated into 
the perturbative quantum chromodynamics. 
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I. Introduction 



The hadronic wave functions in terms of quark and gluon degrees of freedom play an 
important role in the quantum chromodynamics (QCD) predictions for hadronic pro- 
cesses. In the perturbative QCD theory (pQCD) 1 , the hadronic distribution amplitudes 
and structure functions which enter exclusive and inclusive processes via the factoriza- 
tion theorems 2-3 at high momentum transfer can be determined by the hadronic wave 
functions, and therefore they are the underlying links between hadronic phenomena in 
QCD at large distances (non-perturbative) and small distances (perturbative). If the 
hadronic wave functions were accurately known, then we could calculate the hadronic 
distribution amplitudes and structure functions for exclusive and inclusive processes 
in QCD. Conversely, these processes also can provide phenomenological constraints on 
the hadronic distribution amplitudes, the hadronic structure functions, and thereby 
the hadronic wave functions. 

Several important non-perturbative tools have been developed which allow specific 
predictions for the hadronic distribution amplitudes or the hadronic wave functions 
directly from theory. QCD sum-rule technique 4-5 and lattice gauge theory 6-7 provide 
constraints on the moments of the hadronic distribution amplitude. One thus could 
model the hadronic distribution amplitudes by fitting the first few moments in terms of 
Gegenbauer polynomials which are the solutions to the evolution equation of hadronic 
distribution amplitude 2-3 . The Chernyak-Zhitnitsky (CZ) distribution amplitudes 4-5 
(dubbed wave functions in the original papers) constructed in this way are good in re- 
producing the correct sign and magnitude as well as scaling behavior of the pion, proton 
and neutron form factors. However, whether the CZ-like distribution amplitude is the 
correct pion distribution amplitude is still an open problem and its correctness should 
not be judged by only its success in reproducing the correct magnitude of the pion 
form factor. Some earlier lattice Monte Carlo calculations 6 , designed to compute the 
pion distribution amplitude directly, were unable to distinguish between the asymp- 
totic form and the CZ form. In a recent improved lattice QCD calculation 7 the second 
moment of the pion distribution amplitude was found to be small than previous lattice 
calculations 6 and the sum-rule calculations 4-5 , and this suggests that the pion distri- 
bution amplitude may close to the asymptotic form rather than the CZ form. From 
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another point of view, as different wave functions may give a same distribution am- 
plitude, there are still ambiguities about the wave function even if we know the exact 
form of the distribution amplitude. Hence it is still necessary to develop methods which 
specify the hadronic wave functions directly. 

In principle, the Bethe-Salpeter formalism 8 and the more recent discretized light- 
cone quantization approach 9 could determine the hadronic wave functions, but in prac- 
tice there are many difficulties in getting the exact wave functions at present 10-11 . One 
useful way is to use the approximate bound state solution of a hadron in terms of the 
quark model as the starting point for modeling the hadronic valence wave function. 
The Brodsky-Huang-Lepage (BHL) prescription 3 of the hadronic wave function is in 
fact obtained in this way by connecting the equal-time wave function in the rest frame 
and the wave function in the infinite momentum frame, and the corresponding wave 
function gives a distribution amplitude which is significant different from the CZ-form. 
In order to give the CZ-like distribution amplitude which is required to fit the exper- 
imental data, a phenomenological model for the hadronic valence state wave function 
has been proposed 12-14 by adding a factorized function to the BHL wave function. The 
distribution amplitude is almost the same as the CZ distribution amplitude except for 
the end-point regions. Recently a light-cone quark model approach of hadrons 15 has 
received attention for the reason that the model can simultaneously fit low energy 
phenomena, the measured high momentum transfer hadron form factor, and the CZ 
distribution amplitudes. The hadronic wave function in this approach, as will be shown, 
is significant different from the factorized wave functions in Refs. 12-13 though both of 
them give the similar CZ-like distribution amplitudes. 

The purpose of this paper is to discuss the consequence from physical constraints on 
the hadronic wave functions, using pion valence state wave function as an example. In 
Sec. II, we will list several general constraints on the pion valence state wave function 3 ' 16 . 
We then analyse, in Sec. Ill, several existing pion wave functions used in light-cone 
formalism of perturbative QCD. We hope to find a wave function which could give 
both the approximate CZ distribution amplitude required to fit the pion form factor 
data and the reasonable valence state structure function which does not exceed the 
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pion structure function data simultaneously, with several constraints on the pion wave 
functions also satisfied. We find, unfortunately, that the present model wave functions 
have failed the requirements. Furthermore, the most recent light-cone quark model 
wave function given in Ref.15 violates the general constraints severely. For example, it 
gives a probability of finding the valence Fock state in a pion much larger than unity 
if the correct normalization of the pion distribution amplitude is retained. Hence this 
wave function suffers from serious flaws. However we show that the power-law form 
of the pionic wave function can satisfy the requirements and it provides a possible 
example to find a good wave function to be used in QCD theory. In order to properly 
evaluate the effect from Melosh rotation 17-18 connecting the rest frame equal-time wave 
function and the light-front wave function in the light-cone formalism, we re-construct 
in Sec. IV the light-cone quark model of the pion in light-front dynamics. We find that 
the contributions from Melosh rotation in the ordinary helicity (A1 + A2 = 0) component 
wave function by using reasonable parameters seem to have no significant effect on the 
calculated distribution amplitude and structure function in comparison with those from 
the BHL wave function. However, two higher helicity (A1 + A2 = ±1) components which 
were not taken into account in previous perturbative QCD emerge naturally in the full 
light-cone wave function as a consequence from the Wignier rotation. In Sec.V we 
present the summary and some comments. 

II. Constraints on the Valence-state Wave Function 

In this paper, we employ the particularly convenient light-cone formalism 1-3 in 
which the description of the hadronic wave functions is given by a set of n-particle 
momentum space amplitude, 



defined on the free quark and gluon Fock basis at equal " light- cone time" r = t + z 
in physical "light-cone" gauge A + = A + A 3 = 0. Here Xi = kf /p + , with J2i x % — 1> 
is the light-cone momentum fraction of quark or gluon % in the n-particle Fock state; 
k±i, with J2i k±i = 0, is its transverse momentum relative to the total momentum 
and \i is its light-cone helicity. Any hadron state can be expanded in terms of this 




(2.1) 



4 



complete set of Fock states at equal r: 

I H >= E„,A t I[dx}[d 2 k±}^ n (xi, k ±i , A;) 11/ 



u(x i p+,x i p ± +k ±t )x l 
ermions Acj 



n €(a?iP + ,a;iPj_+fc±i)A i i 
gJuons ,/=r I n > 



(2.2) 



with the normalization condition, 

£ | [dx] [d 2 £j | ^(x 4 , fc ± ., Aj) | 2 = 1 , (2.3) 

n,\ 

where the sum is over all Fock states and helicities, and 

n n 

[dx] = 5(1 - ^2xi)]]_ dxi , 

i=i i=i 

n n 

[d 2 k ± ] = 16 7T 3 5 2 (^2 k±i) \\{d 2 A;_Lj/167r 3 ) . 
i=i i=i 



The quark and gluon structure functions G q /h(x,Q) and G 9 /h(x,Q), which control 
hard inclusive reactions, and the hadron distribution amplitudes (/>h(x,Q), which con- 
trol hard exclusive reactions, are simply related to these wave functions 1-3 : 

G a/H {x,Q) = ^ [d 2 k±_i][dxi] | if> n (xi, k±i) | 2 6(x-x a ), a = q or g ; (2.4) 

n J 

and 

rQ 2 - 

(f) H (x,Q) = J [d k± i ]i/) valence (x i , k ±i ) . (2.5) 

In the case of inclusive reactions all of the hadron Fock states generally participate; 
whereas in the case of exclusive reactions perturbative QCD predicts that only the 
lowest particle number (valence) Fock state contributes to the leading order in l/Q. 

In principle the hadronic wave functions determine all properties of hadrons. From 
the relation between the wave functions and measurable quantities we can get some 
constraints on the general properties of the hadronic wave functions. In the pionic 
case two important constraints on the valence state wave function (for the Ai + A2 = 
components) have been derived 3 from n — > /m> and n° — * 77 decay amplitudes: 

£dx J(d 2 k±/I6ir 3 )^(x,k ± ) =U/2y/3 , (2.6) 



and 



C dx rj} qq {x, k± = 0) = , (2.7) 

J 



where f n ss 93 MeV is the pion decay constant. The Ai + A 2 7^ component valence 
state wave functions do not contribute in the tt — > fiu and 7T° — > 77 processes thereby 
the presence of the Ai + A 2 7^ components does not alter the above two constraints. 
Experimentally the average quark transverse momentum of the pion, < k\ > n , is of 
the order (300MeV) 2 approximately 16 . The quark transverse momentum of the valence 
state pion, defined to be 

< A >m= J {d 2 k ± /16n 3 )dx | k\ \\ ^{x, k ± ) \ 2 /p qq , (2.8) 

should be large than < k 2 > n . We thus could require that \J < k 2 > qq - have the value 
of about a few hundreds MeV, serving as the third constraint. The fourth constraint 
is the most natural one: The probability of finding the qq Fock state in a pion should 
be not larger than unity, 

P q q = J (d 2 k ± /16n 3 )dx I ^{x, k ± ) | 2 < 1 . (2.9) 

For the distribution amplitude we compare the calculated one with the CZ-form, 

0cz(^) = 5v / 3 f n x(l-x)(2x-l) 2 . (2.10) 

As to the case of structure function, it should be noticed that the valence state structure 
function is only one part of the valence structure function of a pion. Hence it is 
reasonable to require that the calculated valence state structure function not exceed the 
structure function data for the pion. We will use the NA3 parameterization 19 of the pion 
structure function in comparison with the calculated valence state structure function. 
It should be indicated that the Q 2 corresponding to the NA3 data 19 is very large; i.e., 
Q 2 = 25(GeV/c) 2 . For Q 2 of the order of a few (GeV/c) 2 the structure function should 
increase at large x and decrease at small x in considering the contributions from QCD 
logarithmic evolution 20 , the higher twist effects and other power-law type sources 21 for 
Bjorken scaling violations, with the shape and magnitude not changed too much. 

III. Analysis of Several Existing Wave Functions 

The hadronic wave function depends essentially on the non-perturbative QCD the- 
ory, and it exhibits the full complexity of non-perturbative dynamics. It is necessary 
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for us to use both theoretical tools and phenomenological constraints in studying the 
hadronic wave functions. In order to get clear on where assumptions and approxima- 
tions have been made and where problems may occur, we need to review some previous 
results in the following analyses of several existing pionic wave functions. 

A. The Brodsky-Huang-Lepage Prescription 

Brodsky-Huang-Lepage suggested 3 a connection between the equal-time wave func- 
tion in the rest frame and the light-cone wave function by equating the off-shell prop- 
agator e = M 2 — (£™ =1 h) 2 in the two frames: 

e = f M 2 - (£? =1 qff , £?=i q t = [CM.] 

' ' \ M 2 -Er=iPi + ™ 2 )M] 5 Er=i^ = 0, ££=i a* = 1 [L.C] 

from which one obtains, for two-particle system with m 1 = m 2 (i.e., g° — <?!>); 

- 2 k± + m2 2 In r>\ 

q "MT^y - m • (3 - 2) 

Then for two-particle state there is a possible connection between the rest frame wave 
function ?/>cm (<?*), which controls binding and hadronic spectroscopy, and the light-cone 
wave function i/)lc(x, k±) by 

fe(9 2 )— -m 2 ) . (3.3) 

As an example, the wave function of the harmonic oscillator model in the rest 
frame was obtained from an approximate bound state solution in the quark models for 
mesons 22 

^CM(q 2 ) = A exp(-q 2 /2(3 2 ) . (3.4) 
By using the connection (3.3) one gets the light-cone wave function, 

il>(xi,k ± ) =A ex P [-^{^- + 

The parameters can be adjusted by using the first three constraints, i.e., Eq.(2.6)- 
(2.8), in Sec.II, 

m = 289 MeV; (3 = 385 MeV; A = 0.032 for < %l >^ (35QMeV) 2 . 



Thus we get the pion distribution amplitude, valence state structure function, and the 
probability P q q, 

4>{ X ) =j (d 2 k ± /167T 3 )^(x,k ± ) 

= 2^ X (1 - x ) eX P[~ 8(3 2 T(l-x)] ^ 



F?{x) =Ea= q ^elG a/H (x) 



p q - = J (d 2 k ± /lQir 3 ) J dx | ip(x, ki 

_ A 2 /3 2 -1 /- s r m 2 



- -4i?r- Jo ^(i - a) e^[- 4^) l & ( 3 - 8 ) 
w 0.296 

Fig. 1(a) presents the BHL distribution amplitude in comparison with that of the 
CZ-form and the asymptotic (AS) form 3 , 

(f)As(x) = V3 f w x(l - x) , (3.9) 

which is the leading term of the evolution equation of the pion distribution amplitude 
for sufficient large Q 2 . It can be seen from Fig. 1(a) that the BHL distribution amplitude 
is very close to that of the AS-form while it is significantly different from that of the 
CZ-form. By using this distribution amplitude one accounts for only 50% of the pion 
form factor data. So it is not the "good" distribution amplitude required to fit the 
data. However, we see from Fig. 1(b), where the NA3 parameterization of the pion 
structure function and the calculated valence state structure function are presented, 
that the calculated F}f(x) seems to be a "reasonable" valence state structure function. 
Therefore the BHL wave function cannot give both a reasonable valence state structure 
function and a good distribution amplitude simultaneously. The fact that the BHL 
wave function did not give a "good" distribution amplitude is why the CZ distribution 
amplitude has received attention since it appeared. 

B. The factorized wave functions 

In order to fit the experimental data and to suppress the end- point contributions 
for the applicability of perturbative QCD, a model for the pion valence wave function 
has been proposed in Refs.12-13 by simply adding a factorized function S(x) to the 
BHL wave function, with S(x) specified by 

S(x) = ( Xl -x 2 ) 2 = (l-2x) 2 , (3.10) 
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where x\ — x and x 2 = 1 — x. It leads a distribution amplitude 

«) = l,( 1 -,)( !l - 1 f np| -^ R ] , (3.11) 

which is of the similar shape as that of the CZ-form except for the end-point regions, 
as shown in Fig. 2 (a). The parameters are adjusted by the first three constraints, 

m = 342 MeV; (3 = 455 MeV; A = 0.136 for < kl >^ (343MeV) 2 . 

and the probability P q q is 0.364, which satisfies the fourth constraints, i.e., Eq.(2.9). 
However, the valence state structure function in this case is 

F 2 V (X) = Ea=q,q G a/H (x) . _ 

= IS^a - - !) 2 , 

which is presented in Fig. 2(b). One sees from Fig. 2(b) that there is an unreasonablely 
large hump in the calculated F 2 v (x). Thereby the factorized wave function (3.11) 
though is "good" in giving the CZ-like distribution which fits the data well, it is "bad" 
in giving a reasonable valence state structure function. 

One may specify the factorized function S(x) by other possible distribution 
amplitudes 2 ' 3 which also are constrained by the first few moments given by sum rules 
in terms of Gegenbauer polynomials. The distribution amplitudes in Ref.23 do not 
have the deep dip at x — 1/2, hence they are different from that of the CZ-form. Fig. 3 
presents the calculated distribution amplitudes and valence state structure functions 
by using the factorized wave functions with S(x) specified by 23 

S(x) = 1 + 0.44 Cl /2 {x 1 -x 2 ) + 0.25 C\ ,2 {x 1 - x 2 ) , (3.13) 

and 

S(x) = 1 + (2/3) C 3 2 /2 ( Xl - x 2 ) + 0.43 Cl /2 {x 1 - x 2 ) , (3.14) 

respectively, where C 3//2 (^) is the Gegenbauer polynomials. The parameters are also 
fixed by the first three constraints in Sec. II. We see from Fig.3(a) that the distribu- 
tion amplitudes are more broad than the AS-form. However, the calculated valence 
state structure functions still have unreasonablely large humps, as shown in Fig.3(b). 
Therefore the factorized wave functions also cannot give both reasonable valence state 
structure functions and good distribution amplitudes simultaneously. 
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C. The light-cone quark model wave function 

A light-cone quark model wave function for the pion was given by Dziembowski 
and Mankiewicz (DM) 15 and it has received attention for the reason that it can fit 
the static properties, the form factor, and the CZ distribution amplitude for the pion 
simultaneously. We indicate that this wave function, though is good in "shape", has 
serious problems in "magnitude". 

The main idea in Ref.15 is reasonable: When one transforms from equal-time 
(instant-form) wave function to light-cone wave function, one should consider, besides 
the momentum space wave function transformation such as the Brodsky-Huang-Lepage 
prescription, also the Melosh transformation relating equal-time spin wave functions 
and light-cone spin wave functions. The following assumptions and approximations 
were made 15 : 

1. It was assumed that the ground-state in the pion is described by the harmonic 
oscillator wave function, and adopted the Brodsky- Huang-Lepage prescription 
of the momentum space wave function. 

2. They made a mock- meson assumption that mesons are a collection of quarks with 
all binding turned off and a mock-meson mass equals to the mean total energy 
of the free quarks. 

3. They adopted the mock-meson mass M M = rn^/A + 3m p /4 = 612 MeV rather 
than the real mass m n = 139 MeV or the mean total energy of the two quarks 
M M = 2 < (m 2 + k 2 fl 2 > > 2(m 2 + < kf >) x / 2 = 950 MeV for the pion 
(with < kl >^ (345Me\/) 2 ). 

4. An approximation k° + m ss 2m was made in the obtained light-cone wave func- 
tion, where k° = (m 2 + k 2 ) x l 2 . 

5. It was impliedly assumed that the pion's helicity to be the sum of the light-cone 
helicity of quarks. 

They got, upon the above assumptions and approximations, the A x + A 2 = com- 
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ponent light-cone wave function for the pion, 

^> = ^(f^l-sSSl • (3 ' 15) 

where a± = xMm + m and 02 = (1 — x)Mm + m. From this wave function one obtains 
the distribution amplitude, 

<Kz) = ^—[aM^xix,) - W 2 x lX iY\ expl- 8p2x ™_ x) } ■ (3.16) 

It was argued 15 that the parameters m = 330 MeV and (3 = 450 MeV are rea- 
sonable at large momentum transfer. The distribution amplitude for these parameters 
is very close to the CZ-form, as can be seen from Fig. 4. The parameter A should be 
fixed by the first constraint in Sec. II. Then one can find, from (3.15), the probability 
of finding the qq(Xi + A 2 = 0) Fock state in a pion, 

Pqq = 3.24 

which is unreasonablely large. If we require that P qq - be less than unity, then the 
magnitude of the calculated distribution amplitude will be less than 30% of that of the 
CZ distribution amplitude, and this distribution amplitude surely cannot fit the pion 
form factor data. The calculated structure function is, 

W =1 T6^[ 2 ( 4 ^i^) 3 -2«i« 2 ( 4 ^ 2 ) 2 (317) 
+a 2 l a 2 2 (4(3 2 x 1 x 2 )} exp[- 4f32 ^_ x) ] . 

Fig. 4(b) presents the comparison of the reduced (by a factor 0.1) valence state structure 
function with the NA3 parameterization of the pion structure function. It can be seen 
that the calculated F 2 V is almost two orders of magnitude larger than the experimental 
structure function data as x — > 1. Thereby the pion wave function (3.15) though is 
good in reproducing the shapes of the CZ distribution amplitude and reasonable valence 
state structure function, it has serious problems in producing the correct magnitude. 

The probability may be less than 1 if we change the parameters. For example, we 
find Pqq ~ 0.86 in the case of m — 330 MeV and (3 = 330 MeV, which were used by 
many authors in studying the low momentum properties for hadrons in the constituent 
quark model framework 24 . However, the calculated distribution amplitude, as shown 
in Fig. 4(a), differs significantly from the CZ-like distribution amplitude. 
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D. The power-law wavefunction 25 

It is interesting to look for a wave function that can satisfy the four constraints: 

u 



(1) /; dxi <K Xi ) - 2 ^ 

(2) /' d Xl if>( Xi , k ± = 0) = ^ 



(3) P M - = J d Xl [d 2 k ± ] | ^(Xi,k ± ) | 2 < 1 

(4) = y [d 2 A;_i_] | i>(xi,k±) \ 2 = at or below data. 

Further, we should examine the electromagnetic form factor F W (Q 2 ) and the average 
transverse momentum < k\ >. 

Now let us consider the power-law form of the pionic wave function 

where L and (5 are some constants, £ = x\ — x 2 , and the normalization 
is chosen so that 

For definiteness, let is CZ distribution amplitude, i.e. 

<t> = sVs Ufa -x 2 ) 2 

Then from the constraints one can get the following results: 

9 L - 1 9 
Pas = — < — for L > 1 

99 14 2L - 1 - 28 J 

G f = f fh- I jfe ) |2 = 45^^11^(1 _ x )(2x - l) 4 
u/n J 16tt 3 1 v ; 1 2L - 1 v A ; 

/40~ 

< /c_l 2 > 1/2 = W — Ti f n ~ 356 MeV independent of L 



In order to fit the structure function data at x — > 1, we have a constraint that L is 

59 

58 



below some number. For example, we can put L = || since a simple valence quark 



distribution function like 

G v/w = |x- 1 / 2 (l-x) 
12 



gives tolerable agreement with data for x — > 1. 

The example tell us it is possible to find wave functions that do succeed and these 
are the type that should be used in any future calculation. 

IV. The Revised Light-cone Quark Model Wave Func- 
tion 

One can easily find that the assumptions 2 and 3 in Ref.15 are in fact inconsistent. 
The unreasonable large P q q for the wave function (3.15) should be an indication of 
some unreasonable assumptions and approximations made in Ref.15. Therefore we 
re-construct the light-cone quark model wave function for the pion based upon the 
Kondratyuk- Terent'ev work 18 on the relativistically invariant wave function of two- 
particle system in light-front dynamics. 

We start our discussion from the SU(6) instant-form (T) wave function for the pion 
in the rest frame (q± + q 2 — 0), 

MQ 2 ) = A exp(-q 2 /2f3 2 )( X \x l 2 ~ x[xl)/V2 , (4.1) 

in which x\^ is the two-component Pauli spinor and the two quarks have 4-momenta 
Qi = (<?°> Q) an d <?2 = (q°i with q° = {m 2 + q 2 ) 1 ^ 2 , respectively. The instant-form 
spin states | J,s >t and the front-form (F) spin states | J, A >f are related by a 
Wigner rotation U J (Ref.26), 

| J,A> F =]T U J aX I J,s> T , (4.2) 

s 

and this rotation is called Melosh rotation for spin-1/2 particles. One should transform 
both sides of (4.1) simultaneously to get the light-cone wave function for the pion. For 
the left side, i.e., the pion, the transformation is particularly simple since the Wigner 
rotations are reduced to unity. For the right side, i.e., two spin-1/2 quarks, each particle 
instant-form and front-form spin states are related by the Melosh transformation 17-18 , 

X \T) =w[(q+ + m)xHF)-q R X l (F)] ; 
X l iT) = w[{q+ + m)xKF) + q L x\F)] , 
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(4.5) 



where w = [2q + (q° + m)]~ 1//2 , q R ' L = q 1 ± iq 2 , and q + = q° + q 3 . We also adopt the 
Brodsky- Huang-Lepage prescription for the harmonic oscillator momentum space wave 
function transformation. Then we get the light-cone (or front-form) wave function for 
the pion, 

2 7* 2 

Mx,k ± )=Aexp[- ™ ^_ J E C^x,k ± ,X 1 ,X 2 )x^(F)xHF) , (4.4) 

x{i x) XiM 

where the component coefficients Cq(x, kj_,X±, A 2 ) for J = 1, when expressed in terms 
of the instant-form momentum q^ = (q°, q), have the forms, 

qfM ± ,U) =^ 2 [(g 1 + + m)(g 2 + + m)-g J 2 ]/ v / 2 ; 

C F (x, fc ± , I, T) = -w^Kqt + m)(q£ + m) - q?}/^2 

C^(x, k ± , T, T) = w lW2 [(qi + m)q^ - (q+ + m)g^]/y/2 

C£(x, k ± , I, |) = Wl w 2 [(qt + m)q? - (q+ + m)q*]/s/2 

which satisfy the relation, 

E C ( f(x,^,A 1 ,A 2 )*C F (a;,^,A 1 ,A 2 ) = l . (4.6) 

Ai,A2 

It can be seen that there are two higher helicity (Ai + A 2 = ±1) components in the 
expression of the light- cone wave function for the pion besides the ordinary helicity 
(Ai + A 2 = 0) components. These two higher helicity components arise from Wigner 
rotations (or Melosh rotations in a strict sense) 27 . We also indicate that one should 
express the instant-form momentum q = (q 3 , q±) in terms of the light-cone momentum 
k = (x, k±) in the calculation. However, the relation between q and k is by no means 
unique, thus in practice one needs to construct models relating them. This leads us 
to discuss the following two possible schemes relating the instant-form momentum 
q — (<? 3 ? q±) an d the light-cone momentum k — (x,k±). 

a). Scheme one 

A reasonable connection between q and k was given in Ref.28-29, 

x = (q° + q 3 )/M ; 
k± = q± , 



in which M satisfies 



, k , 2 + m 2 



M = . (4.8) 

x(l-x) V ; 
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From Eq.(4.7) we find, 

q° = [xM + (m 2 + kl)/xM}/2 ; 



(4.9) 

q 3 = [xM - (m 2 + k 2 )/xM}/2 , 



thus 



We notice 



q + = xM , 

- 4.10 
2q + (q° + m) = (xM + mf + kf . 

which is consistent with Eq.(3.2) in the Brodsky-Huang-lepage prescription 3 . The 
detailed reasons for the connection (4.7) can be found in Refs. 28-29. Thus the Ai + A 2 = 
component wave function can be obtained, 

Hx,k ± ) = A T 2 - k K exp[- ™ 2 + k * J , (4.12) 

^ K ' \{a 2 + k 2 ){a 2 + k 2 )YI 2 ^ 8(3 2 x(l-xY 1 ; 

where a\ = xM + m and a 2 — (1 — x)M + m. The constraint P qq < 1 will be satisfied 
with any reasonable m and f3 by using the first constraint to fix the parameter A. 
Fig. 5 presents the calculated distribution amplitudes and structure functions with two 
sets of parameters, i.e., m = 330 MeV, f3 = 540 MeV, and m = (3 = 330 MeV, 
respectively, for the wave function (4.12). When compared with the results from the 
BHL wave function (3.15) in Fig.l, we find, in contrary to the claims in Ref.15, that the 
effect from Melosh rotation seems to have no significant effect on both the calculated 
distribution amplitude and the valence state structure function. The corresponding 
P q q for the two set parameters are 0.228 and 0.552, respectively. 

b). Scheme two 

We adopt the assumptions 1 and 2 in Sec. Ill as the starting point for scheme two. 
In this case the relation between q and k is very simple, 

x = (q° + q 3 )/M ; 
k± = q± , 

where M, the mock-meson mass, defined to be the mean total energy of the free quarks, 
M = 2<(m 2 + g 2 ) 1 / 2 >^2(m 2 + 3/2<£ L 2 >) 1 / 2 , (4.14) 
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the value of which should be, approximately, 1130 MeV and 910 MeV for the above 
two sets of parameters with < k 2 > being (374 MeV) 2 and (256 MeV) 2 respectively 
in scheme one. From Eq.(4.13) it follows, 

q + = xM ; 

2q + (q° + m) = (xM + m) 2 + kl , 

Thereby the Ai + A2 = component wave function should be (4.12) with fixed M 
rather than M in Eq.(4.8). We find P qq to be 0.484 and 0.723 respectively, which are 
larger than those in scheme one. The calculated distribution amplitudes and structure 
functions, as presented in Fig. 6, are also larger than those in scheme one. We also 
calculate the distribution amplitude and the valence state structure function for the 
first set parameters by using M = 612 MeV and compared them with the results in 
Ref.15 as shown in Fig. 4. It can be seen that the distribution amplitude is close to the 
CZ-form, with also the unreasonablely large valence state structure function as that 
in Sec. III. C. The probability P q q is 2.93, which is larger than unity. Thus the flaws 
suffered by the wave function (3.15) mainly raised from the inconsistent assumption 3; 
i.e., the use of a smaller M. The approximation 4 also has a large consequence on the 
results, thereby this approximation seems to be too strong. We also notice that 

q 2 = \[xM + (^±tl)] 2 -m 2 , (4.16) 

which is inconsistent with Eq.(3.2). Further observation of the unreasonableness of the 
results in Ref.15 is that the effect from Melosh rotation should disappear for k± = 0. 
This aspects is satisfied for the wave function (4.12), whereas it is not satisfied for the 
wave function (3.15). 

We only comment without argument here that there are ambiguities of introduc- 
ing a factor, such as ^Jl/2xix 2 adopted in Ref.28 or yj 1 M/AX1X2 adopted in Ref.29, 
to the Brodsky-Huang-Lepage wave function as a consequence from the jacobian re- 
lating instant-form momentum and light-cone momentum. However, the qualitative 
conclusions in above analyses will not changed though the quantitative results will be 
different if the factor is introduced. 



16 



V. Summary and Comments 



After the analysis of several existing pion wave functions in light-cone formalism 
and the re-construction of the light-cone quark model wave function, it is shown that 
the present wave functions (Ai + A2 = component) have failed the requirements to 
fit the pion form factor data and the reasonable valence Fock state structure function 
which does not exceed the pion structure function data with reasonable parameters. 
However, as an example, we consider a power-law form of the pionic wave function. 
It is shown that all of constraints can be satisfied as L — ||. This means that it is 
possible to find wave functions that do succeed and these are the type that should be 
used in any future calculation. 

We also find, in contrary to pervious claims, that the effect from Melosh rotation in 
the ordinary helicity (Ai + A2 = 0) component wave function cannot reproduce a CZ- 
like distribution amplitude. However, as mentioned in Sec. I, it is still an open problem 
whether the CZ-like distribution amplitude is the correct pion distribution amplitude, 
and the large second moment of the CZ distribution amplitude are not reproduced by 
a recent improved lattice QCD calculation. Thus we could consider the quark model 
evaluation of this paper as a suggestion that the pion distribution amplitude may close 
to the asymptotic form rather than to the CZ form. 

Then we meet the problem that the "naive" asymptotic pion distribution amplitude 
can only account for about 50% of the existing pion form factor data. In previous 
perturbative QCD work the adoption of a CZ-like distribution amplitude was found 
to be a possible way to resolve this problem. We speculate in this paper that the 
introduction of the higher helicity states may provide an alternative way to address 
the problem concerning the applicability of perturbative QCD. These higher helicity 
components emerge naturally in the full light-cone wave function as a consequence 
from the Wigner rotation 26 relating spin states in different frames. The existence of 
the higher helicity components in the light-cone wave function for the pion is not in 
contradiction with the requirement of angular momentum conservation, since the full 
light-cone wave function, e.g., (4.4), is an eigenstate of the total spin operator I 2 in 
front-form dynamics 29 . 



17 



The introduction of the higher helicity component wave functions into the light- 
cone formalism may have significant consequence in several problems concerning the 
applicability of perturbative QCD. It is speculated that the perturbative contribu- 
tions from the higher helicity states may provide the other fraction needed to fit the 
pion form factor data besides the perturbative contributions from the ordinary helicity 
states evaluated by using the "naive" Brodsky-Huang-Lepage distribution amplitude. 
However, there are still some difficulties in incorporating the higher helicity component 
wave functions into the conventional perturbative QCD framework since the distribu- 
tion amplitudes for these components vanish from Eq.(2.5). Thereby it is necessary 
to developed a technique which enables us to calculate the perturbative contributions 
from these components. Some progress has been made in this direction and it will be 
given elsewhere. 
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Figure Captions 

Fig.l (a). The normalized distribution amplitude (j)(x) = <p{x)/y/% f w : curves CZ and 
AS are the Chernyak-Zhitnitsky distribution amplitude (Ref.4) and the asymp- 
tonic distribution amplitude (Ref.3); curve BHL is the distribution amplitude 
from the Brodsky-Huang-Lepage wave function (3.5) (Ref.3), respectively, (b). 
The structure function: curve NA3 is the NA3 parameterization of the pion 
structure function; curve BHL is the valence state structure function for the 
Brodsky-Huang-Lepage wave function. 

Fig. 2 Similar as in Fig.l. (a). Curve H is the distribution amplitude (3.11) from 
Huang's factorized ansatz (Refs. 12-13). (b). Curve H is the valence state struc- 
ture function for Huang's factorized wave function. 

Fig. 3 Similar as in Fig.l. Curves BF(I) and BF(II) are the results from the factorized 
wave functions, with S(xi — x 2 ) to be the two sets Braun-Filyanov (Ref.24) like 
distribution amplitudes, for the parameters m = 324 MeV, (3 = 432 MeV, and 
m = 346 MeV, (3 = 461 MeV, respectively. 

Fig. 4 Similar as in Fig.l. Curves DM(1) and DM(2) are the results from the Dziembowski- 

Mankiewicz (Ref.15) wave function (3.15) with the parameters to be m = 330 MeV, (3 = 
450 MeV, and m = 330 MeV, (3 = 330 MeV, respectively. 

Fig. 5 Similar as in Fig.l. Curves a and b are the results from the (Ai + A 2 = 0) com- 
ponent light-cone wave function (4.12), with M in Eq.(4.8), for scheme one with 
the parameters to be m = 330 MeV, (3 = 540 MeV, and m = 330 MeV, (3 = 
330 MeV, respectively. 

Fig. 6 Similar as in Fig.l. Curves a, b and c are the results from the (Ai + A 2 = 0) 
component light-cone wave function (4.12), with fixed M, for scheme two with 
the parameters to be: m = 330 MeV, M — 1130 MeV and (3 = 540 MeV; m = 
330 MeV, M = 910 MeV and (3 = 330 MeV; and m = 330 MeV, M = 612 MeV 
and (3 = 540 MeV. 
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